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Abstract 

The emission of electromagnetic waves from a source within or near 
a non-trivial medium (with or without boundaries, crystalline or amor- 
phous, with inhomogeneities, absorption and so on) is sometimes studied 
using the reciprocity principle. This is a variation of the method of Green's 
functions. If one is only interested in the asymptotic radiation fields the 
generality of these methods may actually be a shortcoming: obtaining 
expressions valid for the uninteresting near fields is not just a wasted ef- 
fort but may be prohibitively difficult. In this work we obtain a modified 
form the reciprocity principle which gives the asymptotic radiation field 
directly. The method may be used to obtain the radiation from a pre- 
scribed source, and also to study scattering problems. To illustrate the 
power of the method we study a few pedagogical examples and then, as a 
more challenging application we tackle two related problems. We calcu- 
late the specular refiection of x rays by a rough surface and by a smoothly 
graded surface taking polarization effects into account. In conventional 
treatments of refiection x rays are treated as scalar waves, polarization 
effects are neglected. This is a good approximation at grazing incidence 
but becomes increasingly questionable for soft x rays and UV at higher 
incidence angles. 

PACs: 61.10.Dp, 61.10.Kw, 03.50.De 

1 Introduction 

The principle of reciprocity can be traced to Helmholtz in the field of acoustics. 
It states that everything else being equal the amplitude of a wave at a point 
A due to a source at point B is equal to the amplitude at B due to a source 
at A. With its extension to electromagnetic waves by Lorentz and later to 
quantum mechanical amplitudes ||] , the applicability to all sorts of fields was 
made manifest. Nowadays the principle is regarded as a symmetry of Green's 
functions when the source point and the field point are reversed. This symmetry 
is actually quite general. As shown in [pi the conditions of time-reversal invari- 
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ance and hermiticity of the Hamiltonian are sufficient to guarantee reciprocity, 
but they are not necessary; in fact, reciprocity holds even in the presence of 
complex absorbing potentials. In the case of electromagnetic waves the only 
requirement is that the material medium be linear and described by symmetric 
permittivity and permeability tensors |^ . This excludes plasmas and ferrite 
media in the presence of magnetic fields. 

In the field of x-ray optics the principle was used by von Laue ||] to explain 
the diffraction patterns generated by sources within the crystal, the so-called 
Kossel lines 0. More recently there has been a widespread recognition that 
these interference patterns contain information not just about intensities but 
also about phases and can be thought of as holographic records from which real 
space images of the location of the internal sources can be reconstructed. Thus, 
under the modern name of 'x-ray holography' there has been a considerable 
revival of interest in this subject [||. 

However, powerful as it is, the usual formulation of the reciprocity principle 
suffers from a rather serious drawback: it refers to the exchange of source and 
field points. As a consequence, a careful application of the principle requires one 
to consider the emission of spherical waves which in crystalline media or even 
in the mere presence of plane boundaries, can be surprisingly difficult (recall 
e.g. studying the radiation by an antenna in the vicinity of the conducting 
surface of the Earth Q or of layered media [0). Furthermore, one is typically 
interested in the asymptotic radiation fields so the relevant exchange should 
involve a source point here with a field point at infinity. 

These technical difficulties have not deterred the users of reciprocity from 
using the principle to make valuable predictions, but a high price has been paid. 
The required asymptotic limits are usually taken verbally and no accounts are 
given of where and how spherical waves are replaced by plane waves. Such 
sleights of hand, because skillfully performed, have not lead to wrong results, 
but intensities are predicted only up to undetermined proportionality factors 
and this excludes applications to classes of problems where absolute intensities 
are needed. Moreover one is left with the uneasy feeling that the validity of 
the predictions is justified mostly on the purely pragmatic grounds that for the 
problem at hand they seem to work which, again, limits applications to problems 
that are already familiar. 

The main goal of this paper (section 2) is to obtain a modified form of the 
reciprocity theorem that gives the asymptotic radiation fields directly and that 
accommodates plane waves and both point and extended sources in a natu- 
ral way. Remarkably the resulting expressions, which include all the relevant 
proportionality factors and yield absolute, not just relative intensities, are very 
simple. 

For many problems the Asymptotic Reciprocity Theorem (ART) obtained 
here represents an improvement not only over the usual form of the reciprocity 
theorem but also over the method of Green's functions. Computing the Green's 
function requires solving a boundary value problems for spherical waves in the 
presence of plane boundaries and/or periodic media; this may well be an in- 
tractably difficult problem. Furthermore, a considerable effort is wasted by first 
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obtaining both near fields and far fields and then discarding the uninteresting 
near fields. The ART is a shortcut that discards the near fields before, rather 
than after they are computed. 

To illustrate the power of the method we consider several applications. The 
first three (section 3) are brief pedagogical examples of increasing complexity. 
First the ART is used to calculate the fields radiated by an arbitrary prescribed 
source in vacuum; next as an application to scattering problems we reproduce 
the kinematical theory of diffraction by crystals. The third example, the radi- 
ation by a current located near a plane dielectric boundary, is straightforward 
when the ART is used but not if other methods are used. One must emphasize 
that what is new in these examples are not the results, but the method; the first 
two are standard textbook material, a special case of the third is treated in |^ . 
As a more involved application of the ART, in section 4 we combine ideas from 
the three previous examples to study two other related scattering problems, the 
specular refiection of polarized x rays by a rough surface and by a continuously 
graded surface. 

The technique of the grazing-incidence reflection of x-rays has received con- 
siderable attention from both the theoretical and the experimental sides 
as a means to obtain structural information about surfaces. The effect of sur- 
face roughness on the reflection is taken into account by multiplying the Fresnel 
reflectivity of an ideal sharp and planar surface by a "static Debye- Waller" fac- 
tor. The problem is to calculate this corrective factor. The calculation has 
been carried out in several different approximations. The Rayleigh or Born ap- 
proximation is satisfactory for rough surfaces with long lateral correlation 
lengths but for x-rays the situations of interest generally involve short lateral 
correlation lengths. Here other approximations such as the distorted-wave Born 
approximation and the Nevot-Croce approximation |l5[ are used. For 
variations and interpolations between these two methods see [Cm, and for a gen- 
eralization to surfaces with non-Gaussian roughness and to graded interfaces of 
arbitrary profile see [0. In these treatments (|14 is an exception) the x rays 
are treated as scalar waves. One expects this approximation to hold at graz- 
ing incidence but at higher incidence angles {e.g., for soft x rays) its validity 
becomes increasingly questionable. Using a modified first Born approximation 
Dietrich and Haase took the vector character of the x rays into account but 
they point out that the validity of their approximation is not in general easy to 
assess and they restrict themselves to studying special interface profiles. 

In section 4 we study this problem using a different approximation; we use 
the ART to develop approximations of the Nevot-Croce type There is, of 
course, a trivial polarization dependence that is already described by Fresnel 
formulas for the reflectivity of the ideal flat step surface. The question we 
address here is whether the "static Debye- Waller" factor shows any additional 
dependence on polarization. The final result is remarkably simple: the "static 
Debye- Waller" factor for the specularly reflected vector waves is the same for 
both polarizations and coincides with that for scalar waves. Finally, some brief 
concluding remarks are collected in section 5. 
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Figure 1: (a) In the usual form of the reciprocity theorem the surface S encloses 
the medium, and all sources, (b) For the asymptotic form of the reciprocity 
theorem the connecting field Ec is a radiation field, its source lies outside the 
surface S — S+ + S . 

2 The reciprocity theorem and its asymptotic 
form 

We wish to calculate the asymptotic radiation fields E and H generated by a 
prescribed current J {t, r) located near or within a linear medium, 

Di = EijEj and Bi = fii^Hj . 

We will assume that the tensors (f) and /i^ (r) are symmetric, but otherwise 
the situation remains quite general, the medium may have an irregular shape, 
or be inhomogeneous, crystalline or amorphous, absorbing, dispersive, etc. 

As in the usual deduction of the reciprocity theorem (see e.g., [Q), we con- 
sider a second set of fields E^ and He, which we will call the "connecting fields", 
generated by a source Jc (see fig.|l|a) . For simplicity we will also assume that all 
fields and sources are monochromatic E = _E(r)e~*'^*, J = J(r)e~'"*, etc. For 
linear media this is not a restriction. 

From Maxwell's equations 

V X E = iKB and V x H = -iKD + — J, (1) 
where K = lo/c, one easily obtains the following identity 

V ■(ExH,~E,xH^^^(E,-J--E-J,y 
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which, on integrating over a large volume V bounded by the surface S, can be 
rewritten as 

E X He- Ec^ Hy ds^— J (^Ec- J- E- j}j dv. (2) 

This expression simplifies if one deals with point sources. For example, consider 
oscillating point dipoles Poe~'"' and PcS"*'^*, located at and fc respectively. 
The current density J is given by J = —itopoS (r — fo) e~"^* and is given by 
an analogous expression. Further simplification is achieved if one assumes that 



the surface S is so remote that the surface integral is negligibly small |18 , then 



Ec {to) -Po^ E (re) • pc . (3) 

This is the usual form of the reciprocity theorem; it says that if we know Ec at 
the location of Po we can calculate E at the location of Pc- This elegant result 
takes us a long way toward a final answer for E, but the remaining problem of 
calculating E^ that is, the calculation of how the spherical wave generated by 
Pc is scattered by the medium, can still be too difficult. 

A more useful version of the theorem can be obtained once one realizes 
that the connecting field is merely a tool that codifies information about the 
influence of the non-trivial medium. Above, the field Ec has been introduced 
by first specifying a source Jc, but clearly this is an unnecessary additional 
complication. In fact, since the most convenient Jc is that which results in 
the simplest Ec it is best to focus attention directly on the field rather than 
its source. Thus we move Jc outside the surface S, to infinity (see fig.^) so 
that throughout the volume V the connecting field Ec is a pure radiation field. 
Furthermore, let the surface S itself be so distant that on S itself both E and 
Ec are vacuum radiation fields. Then 

J i^E x{V X Ec) -EcX{V X E)j ■ ds^ J Ec ■ Jdv. (4) 

At this point it is not yet clear that this form of the reciprocity theorem is 
simpler than eq. J|) but one remarkable feature can already be seen: eq. (||) 
relates the field E at a distant surface S to its source J within a nontrivial 
medium without having to calculate E in the vicinity of J. The "connection" 
between the distant radiation field E and its source J is achieved through the 
much simpler (i.e., hopefully calculable) "connecting" field Ec- 

To bring eq. (||) into a form that is manifestly simpler than (||) the surface 
S is chosen as a cube with edges of length L ^ oo. In fig|l|3 the upper face, 
defined by a constant z coordinate, z = has been singled out as S'-|_, the 
remaining seven faces are denoted S . 

On the upper face S+ we write the field E as a superposition of outgoing 
plane waves of wave vector k satisfying k ■ k = liP' j(? — and k^ > 0, 

^(^)= / -—^2,:5{k-K)E{k)e^-\ (5) 
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where, in a self-explanatory notation, k ■ r = k± ■ fj_ + k^zj^. It is here, by the 
very act of writing E in this form, that the asymptotic limit of discarding near 
fields is being taken. For z << z+ additional terms describing the near fields 
should be included. 

The integral over dkz is most easily done using 



K 

6{k-K) = — 
kz 



S [kz - J - kl) - S Ikz + J - kl 



(6) 



The result is 



where kz = +\/K'^ - k\. 

The choice of the connecting field Ec is dictated purely by convenience. A 
particularly good choice for Ec is the superposition of an incoming plane wave 
of unit amplitude (we use " to denote vectors of unit length) and wave vector 
kci with kc ■ kc — and fccz < 0, plus all the waves scattered by the medium, 

Ec{r^^ece'^^-'' + El{r). (8) 

On S+ , the scattered field E^ is a superposition of outgoing plane waves and is 
given also in a form analogous to eq. (|^), 



with k ■ f — kj_ ■ r± + kzZ^ and fc^ = +y K'^ — k'^. 

Now we are ready to calculate the surface integral on the left hand side of 
eq. (|^). Substituting (^) into (^) the integral over 5*+ separates into two terms, 

one due to the incoming plane wave ec e^^"'^ , and the other due to the scattered 
waves E^ (r) . The first term is 

Ii= f dxdyiz ■ (e X {ike X - x (V x E)\ e*^= ^ (10) 

and substituting ^ its evaluation is straightforward. The integral over dx dy 
yields (27r)^(5(fcj_ + fcc±)- Since k^ = k^ = K^, kcz < and kz > Q this impHes 
that the plane waves superposed in (H) yield a vanishing contribution except 
when k = —kc- Thus, 

^2iKec- E{-kc) (11) 
The contribution of the scattered waves E^ (f) , 

1-2= J dxdyez ■ (^E x {V x EJ - E^ x {V x E)^ , (12) 
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is calculated in a similar way. Substitute (0) and (||) and integrate over dxdy 
to obtain a delta function. This eliminates all Fourier components except those 
with k± = — fc_L. Since — — K"^, and both kz,k^ > this implies k^ = 
k[ > 0. Thus, 



fk^K^. 
{2nf K 



E{k) X (^ik X - E^{k) X [ik X E{k)^ , (13) 



where k = — fc + 2k^ez- Further manipulation using k ■ E^{k ) — k ■ E{k) — 
gives Sz ■[■■■] = 0, so that 

h = 0. (14) 
According to eq. ( |ll| ) and (|^ the only contributions to the surface integral 
over the distant plane 5+ come from products of outgoing with incoming waves. 
Products of two outgoing waves yield vanishing contributions. This result ap- 
plies also to the remaining seven faces of the cube S. Since on each of these 
faces there are only outgoing waves we find that the integral over S makes no 
contribution to the left hand side of (^. Incidentally, this argument completes 
our previously unfinished deduction of the usual form of the reciprocity theo- 
rem, eq. ^ : if both sources J and Jc are internal to the surface S the surface 
integral in eq. (||) vanishes because it only involves products of outgoing waves. 

Substituting ( pT| ) into (^) leads us to the main result of this paper, the 
asymptotic reciprocity theorem, 

e • E{k) = - — / E^- Jdv. (15) 

In words: 

The field E(k) radiated in a direction k with a certain polarization e is 
— 27r/c times the "component" of the source J{r) "along" a connecting field 
Ec{r) with incoming wave vector kc — —k and polarization Cc = e. 

Typically one is interested in the intensity radiated into a solid angle df2; 
since the amplitude E{k) that appears in (|^) and ( p^ ) is not quite the Fourier 
transform of E{f) it may be useful to derive an explicit expression for dW/dfl. 
The total power radiated through the plane 5*+ is given by the flux of the time- 
averaged Foynting vector, ^ Re x B*], 

W = [ d^x^^Re [E X B*] • = f dfl^ (16) 

J OTT J d\l 

Using (0) and d'^k± = k±dk±d(f> — KkzdVL (where 4> is the usual azimuthal angle 
about the z axis) we get 

so that 

In the next section we offer a few illustrative examples of the ART in action. 
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3 Some simple examples 



The ART, eq. (|T5|), holds for an arbitrary hnear medium. In particular it holds if 
the medium is vacuum. Our first trivial example is the radiation by a prescribed 
current in vacuum. Next, to show that the ART can be used to study scattering 
problems we deal with another equally trivial example, the kinematical theory 
of diffraction by crystals. The third example, the radiation by currents located 
near a dielectric boundary, is also straightforward. What is remarkable here is 
the ease with which the results are obtained compared to conventional methods 

3.1 Radiation in vacuum 

In this case the connecting field is just an incoming plane wave, Ec {r) — ic e^^"'^ . 
The ART, eq. (|l|), gives the radiated field with polarization e = Cc as 

e.E{k) = -— [ J{r)-ee-'^-^dv = -—e-J{k), (19) 

so that 

E{k) = —kx{kx J{k)). (20) 
c 

The radiated power, eq. ([l8|), is 

(21) 

as expected. (For radiation by a point dipole just substitute J{k) = —icKp.) 



avv 



Sttc 



kx {kx J{k)) 



3.2 Bragg diffraction 

Consider a crystal described by its dielectric susceptibility x(r) ||l^ which for 
X rays is quite small (typically about 10~^ or less). An incident plane wave 
j7]^giko-r induces a current 

J(f) = -zuPif) = :^x{r)Eoe'^"-\ (22) 

which radiates. The connecting field needed to calculate this radiation is a 
simple incoming plane wave, Ec (r) = Sc e^^"''^ , and the ART, eq. (p^), gives the 
radiated field as 

E{k) = -'-^kx{kxEo)x{k-ko). (23) 
Zc 

The scattered field is proportional to the Fourier transform of the susceptibility 
of the medium; for a periodic medium this is Bragg diffraction. 
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3.3 Radiation in the vicinity of a reflecting surface 

Consider a current (r) located within a uniform medium with dielectric sus- 
ceptibility Xq occupying the region z < (see fig. ^). To calculate the radiation 
in the direction k with polarization e we choose as connecting field an incom- 
ing plane wave with wave vector kc = —k and unit amplitude ic ^ e plus the 
corresponding reflected and transmitted waves, 

4 ir) ^ i "^f t " > " (24) 

I ecte*'=<='■^ for z < ^ ' 

The various wave vectors are given by 

kc = — _ftr cos 9 ix — qe-z = , (25) 

kcr — ~K cos 9 Cx + qe-z , (26) 
kct = ~K cos 9ex~ qiz , (27) 
where K = lu/c, and the normal components q and q are given by 

g = if sin 61 and q ^ K {sin^ 9 + Xo)^^^ ■ (28) 

The amplitudes Scr and Sd of the reflected and transmitted waves are given by 
the Fresnel expressions 

Scr = rsBcr whcrc Ts = z : — , (29) 

q + q Cc-ect 

and 

£ct = tsect where = z- — —■ (30) 

q + q Cc- Cct 

(ccr and Bet are unit vectors describing the polarization of the specular reflected 
and transmitted waves.) 

Then for a source Ji„(r) located within the medium, the ART, eq.(pj|), gives 
the radiated field as 

e • E{k) = / J,„(f) • e^'^'- '- dv . (31) 

C J z<0 

On the other hand, had the source Jout{r) been located outside the dielectric 
medium (z > 0) the corresponding radiated field would be 

e . E{k) ^ - — / J„„,(r) • U e'^^-^^ + e^r e'^'^-'^l dv . (32) 

For an oscillating dipole on the z axis, J{r) = —iujpS (r — ZpC^), eq.(|3ll) and 
(H give 

- ^(t\ _ \ 27riif (p • ee^^'JZp +p . e^^r^ e*''"'p) if Zp > . ^ 

''■^^''^ " \ 27rzifp-e,ti,e-*9^- if Zp < 0. ^'''^^ 
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Figure 2: The connecting field for radiation in the presence of a reflecting 
medium includes reflected and transmitted waves. Here the source Jin is shown 
within the medium (z < 0). 



The power radiated with polarization e, eq.(18), is 



dW 



LStt J e • p 



2 



,2iqzp 



for Zp > , 



(34) 



and 



dW 

dn 



2l Get • p 



2 




e • p 



—tqzp 



for z„ < . 



(35) 



In these two expressions we can recognize the first factor (in square brackets) 
as the power radiated by a dipole in vacuum. The second factor accounts for 
the presence of the dielectric medium. 

4 Specular reflection of polarized x rays 

In this section ideas from the three previous examples are combined to study 
two similar and considerably more involved scattering problems, the specular 
reflection of polarized x rays by a rough surface and by graded interfaces. We 
show that within approximations of the Nevot-Croce type grading and roughness 
affect the specular reflectivity in a manner that is independent of the polarization 
of the incident radiation. 



10 



4.1 Reflection by rough surfaces 

The dielectric susceptibility x(r) that describes the rough surface from which 
we wish to scatter x rays is given by 

/ N / for z > Cix,y) 

= 1 Xo for z<ax,y) ^^'^ 

where the height (^{x, y), is a Gaussian random variable with zero mean, {0 = 0, 
and variance (C^) — (see fig. 

To apply the ART it is convenient to rewrite x (^) 

X (r) = Xs (r) + Sx{r), (37) 
where Xs represents a medium with an ideally flat surface at zq, 

Xs irl = ( ° for z > zo 

^ ^ \ Xo for 2: < zo ^ ^ 

and dx represents the roughness. 

Let Co e*'^"''' be the incident field. The total scattered field e(f) includes the 
wave Ss (r ) specularly reflected by the step Xs (r) plus waves de (r) scattered by 

Sx (r) 

sif)^esir^+Se{r). (39) 
The first term on the right is 

e-^f) = e,r,e-2'«^°e'"^'-^ (40) 

where 

q-qir-et , , 

Ts = — — (41) 

q + q eo ■ et 

(Br and it are unit vectors describing the polarization of the specular reflected 
and transmitted waves). The second contribution in eq.(|39|), the field Se{r) 
includes a specular component plus diffusely scattered and evanescent waves, 

6e (f) = Ser e.e'*^"-"^ + 6ed{r). (42) 
Using eq.(0) this may be written as 

Se{r)= /^^fc-(fc)e^^"'n (43) 

where 

6e{k ) = ^<5e. {2nf 5{k^ - fcoi) + 5ed{k)- (44) 

To calculate 5e{r) we can proceed exactly as in the previous section (3.3): 5e{r) 
is the field radiated by a current 5J{r) in the presence of the medium Xs (^)- 
The current 

6J{^ = ^&x{rlE{r), (45) 
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Figure 3: The problem of scattering by a rough surface can be tackled using the 
ART by adding a fictitious overlayer Sx- 



originates in the polarization of the roughness 5x{r) by the total electric field 
E{r) due to the incident and all scattered waves, including those generated by 
the roughness itself. Thus, the challenge here is that the field E (f) is itself 
unknown; an approximation for it must be obtained as part of our solution. 

We can exploit the arbitrariness in the separation of x into Xs phis 
6x (r) to suggest a self-consistent approximation for E. Suppose we choose zq 
positive and considerably larger than the roughness a (see fig. ^). Then Sx 
represents a fictitious overlayer that extends well into the vacuum; the sign of 
6x (r) is opposite to that of Xs ^^^d in the vicinity of zq they completely cancel 
out. The field 6e{k) in a direction k with polarization e is given by eq.(^l[) 

e • de{k) ^'-^ fdv Sxir^E (f) • e^t e^^-'^" , (46) 



where the connecting field is precisely as in eqs.(p4|)-(30) except for phase shifts 
due to the reflecting surface being at zo, 



e 



-2*«^"r,ee. and e^t = e^^'-'^'n.eet. (47) 



The reason behind the somewhat surprising choice for zq will now become 
clear: slightly above zq, in vacuum, the exact field is 

E{r) = eoe'^^-^ + e{r) = eoe'^°-^ + £,(r) + Se{r), (48) 

but slightly below zq and, in fact, over all of the extension occupied by Sx (r), 
we are also in vacuum {6x {r) and Xs cancel each other) and therefore E (r) is 
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given by the same expression (p8|). The last term Se(f), given by (p2|), includes 
some weak diffusely scattered and evanescent waves Sed{r)- Our approximation 
consists of neglecting them. Therefore, 

E (f) « eo e'''"-^ + rere'''--^, (49) 

where the specular reflections by Xs ^^'^ have been combined into the 

single, and still unknown, reflection coefficient r, 

r = r,e-2»9^° + fc^. (50) 

Substituting into eq.(|4^) yields 



e • Se{k) 



dx dy 



dz 



£ct e 



ikct -r 



(51) 



From now on we focus our attention on the specularly reflected component; let 
e ^ Be — er, icr — eo, k — kj. ^ —fee- Substituting eq.(Q) into the left hand 
side {l.h.s.), using {2tt)'^ S{k± — ko±) = (27r)^(5(0) — J dx dy we get 



l.h.s. ^ ^6er {2ttY (5(0) = ^{r- r,e-^"^^°) J dx dy. (52) 

This shows that the unknown reflection coefficient r we want to calculate appears 
in both the left and the right hand sides of (|5^), as part of the radiated field and 
also as part of the field that induces the source; eq.(|l]) permits a self-consistent 
calculation of r. 

The integral over dz in the right hand side of eq.(pT|) is elementary and the 
remaining integral over dx dy is performed using the identity 



/ dxdyCQ^'^''-y^ 
J dx dy 



(53) 



where ^ is a Gaussian random variable with zero mean, (C) — 0, and variance 
(C^) = cr^. The right hand side (r.h.s.) of eq.(^l|) becomes 



r.h.s. 



dx dy 



eo ■ get 
q + q 



q 



gl(<j-q)zo _ g- 



(54) 



which can be further rewritten by substituting Sct as given by eq.(^7|), and using 



r =K^Xo^ and 



e • Cct 



e • Cct 



eo ■ et 



and = 1. 

e • Cct 



(55) 
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Finally, equating eq.(^2|) to ( p^ yields a self-consistent approximation to r. 



g -I- g eo • et 

This coincides exactly with the Nevot-Croce result for the polarization eo = e* = 
Cr for which the ratio eo-ei/ Cr-et is unity, and provides the correct generalization 
to all polarizations. According to this approximation the specular reflection 
coefficient r has no polarization dependence beyond that already implicit in the 
reflection coefficient for the ideal flat step surface; the "static Debye- Waller" 
factor exp{—2qqa^) is polarization independent. 

Notice that any possible dependence on the arbitrary choice of zq has can- 
celled out. 



4.2 Reflection by smoothly graded surfaces 

The problem of scattering by a smoothly graded interface is similar and some- 
what simpler. Here the susceptibility xi^) depends only on the normal coor- 
dinate z and not on the transverse coordinates x and y. This implies that the 
tangential component of momentum is conserved in the scattering; there are no 
diffuse waves, there is only specular scattering. 
As before, it is convenient to separate x (z) into 

x{z)^Xsiz) + Sxiz), (57) 
where Xs (^) represents an ideally flat surface at zq, 

Xs = " ^ "° (58) 

-^"^ ' \ Xo for z < zo ^ ^ 

and dx (z) is an overlayer (see fig]4|) describing the smooth transition from bulk 
to vacuum. ^ 

Let eo e^''"'^ be the incident field. The total scattered field e'(r), eq.(|39|), 

e{r)^Ss{r^+6eir}. (59) 

includes the wave es{'r) reflected by the step Xs (z), eq. (pO|) , plus waves (5e (r) 
scattered by the overlayer 6x (z)- While diffusely scattered waves are not present 
in Se{r), faint evanescent waves could be; these are weak near field effects and 
we neglect them. Thus 

Se{r) =6erere'''^-^, (60) 

and the Fourier expansion, eq. (|4^) , and transform 5e{k), eq. (^4[) , remain other- 
wise unchanged. Once again, Se (f) is radiated by a current 

5M^^Sx{z)E{r), (61) 

where the field E{r) includes the incident and the unknown reflected waves; 
E (r) must be self-consistently obtained as part of the solution. Then the ART, 
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Figure 4: The problem of reflection by a graded surface can be tackled using 
the ART by adding a flctitious overlayer 5x- The hatched region shows the 
transition region from (5^ = to (5x = ~Xo- 

in the form of eq.(|3^), gives the field (5e(fc) in a direction k with polarization e 
as 

e • 5e{k) = '-§Jdv Sx{z)E (r) ■ e^t e^'^- *' , (62) 

with the same connecting field given back in eq.(^). 

The approximation we use for E (r) is the same as in last section. The 
arbitrariness of zq can be exploited by choosing it large enough that the overlayer 
extends well into the vacuum. Near zq the overlayer and the sharp step Xs 
cancel each other out; slightly above zq, in vacuum, the field is 

E (r) = eo e''^"-^ + re^e*'^'-^ (63) 
where r is the unknown reflection coefficient we want to calculate, 

r = + fc^. (64) 

Slightly below zq and over most of the extension occupied by Sx (z) we are also 
in vacuum (provided the bulk to vacuum transition is not too gradual) and we 
approximate -E (r^ by the same expression, eq.(|63|). Substituting into eq.(|62|) 
yields an equation for r, 

^{r~ r,e-2^«^«) (2^)' S{k^ - k„^) = 
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iK 

T 



— CxD 



dz5x{z) dxdy eoe''=«"" + re^e"^'--'' ■ Sd e'""'-'' . (65) 



The integral over dxdy yields a delta function, {2tt) S{k± — kQ±), and we can 
substitute e = Cc — Sr, Scr — eo, k = kr = — fee- The integral over z is 
conveniently expressed as 
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dz Sx(z) e-'^' = ^ [ e-*^"" + ^^^] , (66) 
«Q Xo 

where x'iQ) is the Fourier transform of dx{z)/dz, 

X'(Q) = f^dz^e-^Q^' (67) 

Eq.(p6|) is proved by integrating the left hand side by parts, using Sx{zq) « —Xqj 
and d5x{z)/dz = dx{z)/dz. Using (f — = K^Xo ^^'^ the identities in eq.(p5[) 
the final result is 

r = rs-rrz r- (68) 

Notice that any possible dependence on the arbitrary choice of zq has cancelled 
out. 

This coincides exactly with the scalar wave result and provides the 
correct generalization to all polarizations. Within these approximations the 
specular reflection coefficient r has no polarization dependence beyond that 
already implicit in the reflection coefficient for the ideal flat step surface; the 
"static Debye- Waller" factor is polarization independent. 

To conclude we mention some illustrative examples: 
(a) The error-function profile 



gives 

x'iq + q) _^^2g^a^^ (70) 



x'iq-i) 

the same factor obtained in the previous section for a Gaussian rough surface. 
This is as expected, the error function is the averaged profile for the Gaussian 
rough surface. 

(b) The Epstein (or Fermi distribution) profile ||l^ 
gives 

x'iQ + q) ^ g + g sinh[7rcr(g - q)] 
X'(g-g) q-q smh[TT a {q + q)Y 
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(c) The triangular profile 




gives 

x'iQ + <?) _q~q sm[{q + q)(T/2] 



(73) 



(74) 



X'iq-q) q + q sin[{q- q)a/2]' 
The reliability of these approximations was studied in in the case of 
scalar waves. There is no reason to expect any difference from the conclusions 
reached there: the "static Debye- Waller" in eq.(|6^) provides a remarkably good 
approximation for the intensities reflected by interfaces of arbitrary grading 
profile even for transition regions that are quite wide (cr as large as several 
nanometers). The phase of the reflected waves is however more sensitive; eq.(|68|) 
provides a good approximation for more abrupt transitions {a of the order of 1 
nm or less). 



5 Conclusion 



The main result of this work, eq.(15), is an asymptotic form of the reciprocity 
theorem which can be used as the basis for a practical method for calcula- 
tions. The theorem states that the field radiated in the presence of a nontriv- 
ial medium, in a certain direction and with a given polarization, is a suitable 
'component' of the radiating source. This 'component' is to be extracted by 
introducing an auxiliary 'connecting' field which contains the necessary infor- 
mation about the medium. The practical advantage of the method lies in the 
simplifications achieved by systematically avoiding unnecessary calculations; it 
thereby allows one to tackle problems of increasing complexity. 

In forthcoming papers we will further explore the application of the ART to 
the study of the dynamical diffraction of radiation generated by sources within 
a crystal, the so-called Kossel lines. Even this well explored topic has not been 
exhausted. Of particular interest are situations where the Bragg angle lies close 
to 7r/2 and the Kossel cones degenerate into single beams and situations 
where the source location is revealed by the oscillatory 'Pendellosung' structure 
of the diffraction pattern |2ll . Other applications will include a new approach 
to thermal diffuse scattering under conditions of dynamical diffraction ||2^ . 

Acknowledgments. I am indebted to P. Zambianchi and E. Sutter for 
valuable discussions. 
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